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Abstract 

This paper gives a simple way of checking the stabil- 
ity with respect to an arbitrary stability region of a 
family of polynomials containing a vector of parame- 
ters varying within prescribed intervals. It is assumed 
that the parameters appear affine multilinearly in the 
characteristic polynomial coefficients. The condition 
proposed here is simply to check the phase difference 
of the vertex polynomials. This test based on the 
mapping theorem significantly reduces computational 
complexity. We omit mathematical proofs in this 
version. The results can be used to determine var- 
ious stability margins of control systems containing 
interconnected interval subsystems. These include 
the gain, phase, time-delay, H°° and nonlinear sec- 
tor bounded stability margins of multilinear interval 
systems. 

1. INTRODUCTION 

The problem of determining whether a control system 
remains stable or not when some of its parameters 
undergo perturbations within prescribed intervals is 
one of the central problems in robust parametric sta- 
bility. The case where the parameters appear affine 
linearly in the coefficients of the characteristic poly- 
nomial has been solved effectively by the Edge The- 
orem [1] and the Box Theorem [2,3] The case where 
the coefficients are affine multilinear functions of the 
parameters is of current interest (see [4] - [9]). This 
case arises in control systems where transfer function 
coefficients perturb or state space parameters perturb 
or matrix fraction factors perturb. The present paper 
deals with the multilinear case also and shows that a 
sufficient condition for the stability of such systems 
is that the phase difference between various vertex 
polynomials evaluated along the stability boundary 
be less than 180°. This is an “image set” approach 
using the convex hull property of multilinear image 
set [10,11,12]. We show how this result can be used 
in conjunction with the multilinear version of the Box 
Theorem [5] to develop a highly efficient phase test 
for stability where only the Kharito nov vertices ar e 


tested. These conditions which are sufficient condi- 
tions can be tightened by increasing the number of 
vertices and this eventually leads to necessary and 
sufficient conditions. Using these results we show 
how various worst case stability margins such as gain, 
phase, time-delay, H°° and sector bounded nonlinear 
stability margins for multilinear interval systems can 
be found. 

2. NOTATION AND MAIN RESULTS 

Let p = [Pi»P 2 > • • • } pi] denote a vector of real param- 
eters. Consider the polynomial 

5(s, p) := $o(p)+$i(p) a +^2(p) s2 H f^n(p)-5 n * (1) 

wherein the coefficients 6*(p) are affine multilinear 
functions of p, i — 0,1 The vector p lies in 

an uncertainty set 

n:={p | Pi < Pi < P+, *=1.2, •••,*}• (2) 

The corresponding set of polynomials is denoted by 



A := {£(a,p) | P E II}. 

(3) 

Let V denote the vertices of II, i.e., 


V := {p 

| Pi = Pt or pi = Pi , » = 1,2,- 

”,0 (4) 

and 

Ay := {6(3, p) | P eV}. 

(5) 


denotes the set of vertex polynomials . 


Fixing s — s*, we let A (s*) denote the set of points 
£(s*,p) in the complex plane obtained by letting p 
range over II: 

A(a*):={tf(5*,p) I pen}. (6) 

Likewise 

Ay(0 :={*(»•, p) I pev}. (7) 

The convex hull of a set of points V in the complex 
plane is co V . Denote the set of convex combinations 
of the vertex polynomials by 

E := {>M<( a ) + (1 - \i) 6 i 00 I *0 6 [0, 1] and 

£ Ay}. (8) 


(NASA-CR— 192893) STABILITY MARGINS N93-25226 

FOR MULTILINEAR INTERVAL SYSTEMS BY 
WAY OF PHASE CONDITIONS: A UNIFIEO 
APPROACH (Tennessee State Univ.) 

5 p 


Unci as 


Let S C <E denote an open stability region and dS 
the boundary of S. Let A(A) denote the image set 
evaluated at s = A of the set of polynomials A. We 
can now state a first result based on the Mapping 
Theorem of [13]. 

Theorem 1. The set of polynomials A is stable 
with respect to S if if has at least one stable poly- 
nomial and i) 0 ^ A(A) for some A E dS ii) the set 
of polynomials E is stable with respect to S . 

The problem of checking the stability of the set E can 
be solved computationally relatively easily in view if 
the fact that E(A) consists of straight lines whose end 
points are contained in the set Ay(A). To formalize 
this let 6(a) be a polynomial and s = A a point in the 
complex plane. Let ^(A) denote the argument of the 
complex number 6(A); 

5(A) := |5(A)|e>'*‘W. (9) 

For a given set of polynomials T let 

$t( a ) : = bu p I- ( 10 ) 

We can now state the following useful corollary. 


Corollary 1. (Theorem 1) The set of polynomials 
A is stable with respect to S if it contains one stable 
polynomial: 



M( -) 


S* * 

m 








Figure 1. Unity Feedback System 


3.1. Gain Margin 


Consider the family of feedback systems as above with 


M(s, p,q) 


p) 

D{s } q) 


( 11 ) 


where N(s, p) and D(s, q) are multilinear functions 
of p and q and (p,q) varies in a polytope P with 
vertex set Py. The upper gain margin of the family 
of the family of feedback systems with M is defined 
as the largest value L* > 0 of L so that 


<0 + (1 + L)N (s, p) (12) 

remains Hurwitz for all L E [0, L *) and for all (p, q) E 
P. The lower gain margin of the family may be sim- 
ilarly defined by replacing (1 + L) by (1 - L ). Let us 
first define the following sets corresponding to each 
fixed real value of L\ 


1l(u>) q) + (1 + L)N(ju } p)| 

(p>q)eP}. ( 13 ) 


a) 0 ^ A(A) 



boundary of the region S. 

The above result states that A is stable if the net 
phase difference of the vertex polynomials evaluated 
along the stability boundary is less than 180°. We 
shall refer to this as the Phase Condition. The 
Phase Condition can be easily verified by the Seg- 
ment Lemma [13] without sweeping over frequency. 

3. CALCULATION OF STABILITY MAR- 
GINS 

Exact calculation of stability margins for a multilin- 
ear interval control system is computationally diffi- 
cult. In this section, we give simple techniques to “lin- 
earize” the problem so that the Phase Condition can 
be used to estimate various stability margins without 
excessive calculations. 


Let 

2l(w) := co {L>(jt*/,q) + (1 + L)N{ju)^ p)| 


(p.q) e V v }- 

(14) 

Now let 

L(ut) := inf{L | 0 E2l(u;)} 

(15) 

and similarly, 


L(u>) := inf{L \0eX L {u>)}. 

(16) 

The corresponding gain margins are; 


L * inf L{w) 

w 

(17) 

L* := infL(u>) 

u 

(18) 

Theorem 2. The gain margin L * 

of M(s ) is 


bounded from below by L * . 

The point here is that lower bound L * can of course 
be found from the Phase condition, since it involves 
only the vertices. 



3.2. Phase Margin 


and define 


Similarly, the phase margin of the family with M(s) 
is defined as the largest value 0 * > 0 of 9 so that 

D{3,q) + e ie N{a,p) (19) 

is Hurwitz for all 6 € [0,0*) and for all (p,q) G V • 
Consequently, we define the sets corresponding to 


each fixed value of 8: 

X*(w) := {£*(jw,q) + e , ' ff lV(iw,p)| 

(p,q)<E7>} (20) 

I*(w) := co {D(jw, q) + e> e N{ju>, p)| 

(p.q)eTV}. (21) 

And let 

4>{w) := inf{tf | 0el*(w)} (22) 

<j>{J) inf{0 I oe T*(o;)}. (23) 

The corresponding phase margins are: 

</>* := inf</>(u;) (24) 

u 

0* inf<£(u>) (25) 

UJ 


Theorem 3. The phase margin <f>* of M(s) is 
bounded from below by <f >* . 

Again this lower bound can be found using the Phase 
Condition . 

3.3. Time- Delay Margin 



Figure 2. Time-Delay System 


Consider the system with the time delay T. The time 
delay margin of the system is the maximum allow- 
able time delay preserving closed loop stability. It 
can be calculated as before from the convex hull of 
the corresponding image sets. Let us define the sets 
corresponding to each fixed real value of T: 

I T (a») ~ {D{jw,q) + e~ juT N{jw,p) \ 

(p.q)eP} (26) 

ZtM := {D{ju>,q)+e-’ uT N{ju,p) \ 

(p,q) G Tv} 


T{u>) := inf{T | OGXr(w)} (28) 

f(w) := inf{T | 0ei T (w)}. (29) 

The corresponding maximum time delays are: 

T — inf T(u>) (30) 

T := inf T(u;) (31) 

Ui 


Theorem 4. The time delay margin T* of M(s) is 
bounded from below by T * . 

As before T* can be found from the Phase Condition . 

3.4. H°° Stability Margin 

Consider the following configuration with Q(s) fixed 
^(^PiQ) as before and AM representing an un- 
structured perturbation, Let 



then the stability margin of the system is defined 
as follows: 


8U P(p,q)e^ ||^( 5 )[1 + Q( 5 )x>(*!qj] 1 ^ 
1 

Q 3 (*)[n(-,^+Qi(*)JV(*Tp)I ^ 

From the Lemma in [14], we have 

V/3>0, ||(3(a)[l + Af(j,p,q)<3(«)] _1 ||oo 

for all Af(s,p,q) iff 

Q 2 (a)D(s,q) + Qi(a)N(3, p) + / 3e? e Qi{s)D(a , q) 
= [Qa(«) + ^ 6 QM]D{3, q) + Qi(s)N{3, p) 


(27) 





is Hurwitz for all 6 6 [0,2ir), and (p,q) E V ■ Now 
let us define the sets corresponding to each fixed real 
value of a and 9: 

lp,e{u) := {[Q2(ju)+0e’ e Qi{ju>)]D{jw,q) + 

Qi(ju>)tf(jw, p) | (p, q) £ V} (33) 
ip,e{w) co {[t? 2 (jw) + Pe je Qi(ju>)]D(ju, q) + 

Qi(ju>)N(ju, p) | (p> q) E V v }• (34) 

If we let 

;= inf {ft | 0 E Ip % e{ u ) for some ( 35 ) 
(3(u>) inf{/3 | 0 E for some 0} (36) 

The corresponding stability margins are: 

ft* — inf ft(w) and ft* =in(ft(w). (37) 

' ' u> 


Theorem 5. The H 00 stability margin ft* of the 
configuration in Figure 3 is bounded from below by 

F- 

The bound ft* can be found from the Phase Condi- 
tion. 


3.5. Nonlinear Sector Bounded Margin 


+ 


•9- 


M(.) 


Figure 4. Lur6 Problem 


Now let us consider Lure problem. According to [15], 
if <7 (a) is a stable transfer function and 4> belongs to 
the sector [0,Jb] f then a sufficient condition for abso- 
lute stability is 


Re {i +ff ( ja ')} > 0> v w61R - ( 38 ) 

In this problem we want to find the maximum value 
lb* of k so that 

Reji + M(jw,p, q)] > 0, V u£E (39) 


for all M(s, p,q). 


From the theorem given in [15], we know that g(s) = 
is SPR iff the following three conditions are sat- 
isfied: 


1) Re{g(0)} > 0 

2) n(s) is Hurwitz stable 

3) d(s) + jan(s) is Hurwitz stable for all a 6 ft. 


Consider 


1 

k 


+ Af(jw,p,q) 


1 , ^(J.P) 
k D(s,q) 

£>(a,q) + fcJY(s,p) , . 

kD{s,q) 


Suppose the transfer function satisfies the above two 
conditions: 


1) Re{ Pt0 ’^o,qi°' P i } > 0 for a11 k e Ml and 

for (p, q) 6 V 

2) D[s, q) + kN{a,p) € H for all k E [0, k *] and for 

(p.q) ev. 


Let us define the following sets corresponding to each 
fixed value of k and a: 

I t ,„H := {kD(ju,q)+ja[D{ju>,q) 

+kN{ju, p)] | (p, q) EV} 

= { k[D(ju , q) -(- jaN(ju), p)] 

+jaD(ju,q) | (p,q) EV} (41) 
21t,„(w) := co {k[D(ju>, q) + jaN(ju>, p)] 

+jaD(jw, q) | (p, q) E 7V}(42) 

If we let 

Jfe(u>) := inf{fc | 0 E Xjk^u/), for some a} (43) 
fc(a^) := inf{fc | 0 E 2*,«(w)i for some a}, (44) 

the corresponding nonlinear sector bounded margins 
are: 

ib* = inf k(w) and £*=inf£(u;). (45) 

u V 


Theorem 0. The nonlinear sector bounded margin 
k* of M(s) is bounded from below by k * . 

4. CONCLUDING REMARKS 

We have given a simple computational method to 
check the stability and compute various stability mar- 
gins of polynomial families containing parameters 
which appear affine multilinearly in the character- 
istic polynomial coefficients and which vary in pre- 
scribed intervals. Furthermore, the proposed method 
completely eliminates frequency sweeping. Although 
the test is only a sufficient condition we expect it to 
be useful especially in view of the practical necessity 



of maintaining adequate margins of stability. Addi- 
tionally the bounds obtained on the stability margins 
can be made as exact as desired by introducing addi- 
tional vertices. Finally it is worth emphasizing that 
the computations given here avoid construction of im- 
age sets, a task which is enormously unwieldy from a 
computational standpoint. 
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